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Abstract 

We present a method to enhance the ripple structure of the scattered electromagnetic field in the visible 
range through the use of Laguerre-Gaussian beams. The position of these enhanced ripples as well as their 
linewidths can be controlled using different optical beams and sizes of the spheres. 
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INTRODUCTION 



In 1 890, Ludvig Lorenz obtained one of the few analytical solutions of the Maxwell equations 
in inhomogeneous media, the scattering of a plane wave from a dielectric sphere Jl]. Later, Gustav 
Mie rediscovered this result and applied it to metallic spheres [2]. His results explained the change 
in colors of colloidal suspensions of gold nanoparticles of different sizes, also showing for the 
first time the size dependence of localized plasmon resonances in nanostructures. The resulting 
theory is usually called Lorenz-Mie theory (or Mie Theory) and solves the scattering problem of 
an incident plane wave propagating in a homogeneous isotropic medium hitting a homogeneous 
isotropic sphere. However, with the advent of modern computation techniques, this theory has 
been developed further to the Generalized Lorenz-Mie theory (GLMT), which solves the scattering 
problem for any incident electromagnetic field yfl. The GLMT has found application in many and 
diverse fields. We refer the reader to Gouesbet et al. [4] for an extensive review of all the possible 
applications. In particular, it is important to highlight the remarkable impact that this theory has 
had in the development of nanopho tonics. Some very recent examples are the study of photonic 
jets [5], the characterization of coherent perfect absorbers Jg], the prediction of the optical pulling 
force OD or the control of localized surface plasmons [8], or the excitation of multipolar resonances 



However, one aspect of the GLMT that is usually overlooked is its intimate relation to the 
angular momentum (AM) of electromagnetic waves Jioll . The GLMT is naturally described in 

(x) 

spherical coordinates, since the scatterer is a sphere. The multipolar modes A -^ are the solutions 
of Maxwell equations in spherical coordinates where j = 1, 2, .., oo stands for the angular mo- 
mentum of the mode, \m\ < j is the third component of the angular momentum and x = e, m 
accounts for the parity of the beam, m being the magnetic multipole and e the electric one. In- 
deed, it can be proved that J 2 A£ = j(j + 1)A£, J,A$ = mA$, FA#? = (-)'A£J and 
PA^ = (— Y +1 Aj^. That is, the multipolar modes are eigenstates of the AM squared operator 
J 2 , the z component of the AM operator J z , and the parity operator P. The field of AM of light 
was significantly developed when in 1992 Allen and co-workers showed that in the paraxial ap- 
proximation one could find a set of modes which had a well defined value of the third component of 
the AM 0, J z . A particular set of paraxial modes with this property are the Laguerre-Gaussian 
modes (LG q j, with q the radial index and I the azimuthal one). This important finding opened up 
a whole new field allowing for innumerable applications in quantum optics, microscopy, biology, 



optical trapping and astrophysics, just to mention a few [12]. Finally, although the total AM of a 
vectorial field is composed of an orbital part (L) and a spin part (S), i.e. J = L + S, these two 
components cannot physically be separated. That is, they can be split mathematically, but they 
do not correspond to any physical observable. When either of these operators are applied to a 
Maxwell field, the result from that operation is not a Maxwell field 11131 . p501 11141 §16]. 

In this article we explore the use of LG modes to excite spherical dielectric nanoparticles. 
Exploiting the properties of angular momentum and using some analytical techniques we manage 
to enhance the ripple structure of an arbitrary sphere in a lossless, homogeneous and isotropic 
medium. Furthermore, the results presented in this paper will allow for a better understanding of 



recent scattering experiments of silica spheres with Laguerre-Gaussian beams ft 1 311 . 



GENERALIZED LORENZ-MIE THEORY 



In this section we solve the following scattering problem: a monochromatic LG g j beam propa- 
gating in a homogeneous, isotropic and lossless medium impinges on a dielectric sphere made of 
a homogeneous and isotropic material. Since all the material properties can be reduced to a pair of 
scalar numbers (the electrical permittivity e and the magnetic permeability /i) for both the sphere 
and the surrounding medium, the electomagnetic field can be obtained in two steps. First, one of 
the two fields (E or H) is found by solving the monochromatic wave equation. Then, the other 
field is obtained by applying the curl to the first obtained field: 

V 2 E + A; 2 E = (1) 
V x E = ikfiB. (2) 

where k = nu/c with n = y/efl the index of refraction of the medium. The temporal dependence 
of the fields is supposed to be exp(— iwt) in the whole text. The multipolar fields j A^, j 
are precisely a basis of solution for the two equations above |l6|l . Furthermore, as mentioned 
in the introduction, they are specially well suited for problems in spherical coordinates as they 
are rotationally symmetric. Thus, we will use them as a basis to decompose all the fields in the 
problem. The fields we have to consider are the incident (E\ H l ), the scattering (E sca , H sca ) 
and the interior field (E l , H'). The incident field is given (it is a LG p z in this case), therefore its 
decomposition in multipoles can be computed. In general, this decomposition cannot be evaluated 
analytically. That is, if the incident field has a decomposition of the form E* = V . m g^A^ + 
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g^A^l then both and are obtained by computing two cumbersome double integrals 
9jm = Jo Jo™ ^ ' A*im s i n ®dQd(j). However, once this has been done, it has been proved that the 
scattered and the interior fields are readily determined by applying Maxwell boundary conditions 
on the sphere |3|]. Thus, the only technicality of the problem is finding the decomposition of the 



incident LG^ beam into multipoles. This decomposition was found in |17|1 : 



j=\i+ P \ 



(3) 



C jlp = r S m6di l+p)p (e)f q>l {ksm9)de 

•J 



where A is the vector potential associated with the incoming beam, p is the polarization of the 
incoming beam whose value is either -1 for left circular polarization or 1 for right circular polar- 
ization, and where m = I + p stands for z component the total angular momentum. The reduced 
rotation matrix d^ (9) can be found in Rial , and the function f q j(k r ) is related to the Fourier 
transform of the incident field, with k r being the transverse momentum, i.e. k r = ^k 2 + k'*. Note 
that the function Cji p is still defined as an integral, although it is no longer a double one. The 
reason why this happens is that the beams we have chosen as incident field have a specified value 
for J z . Actually, eq. © is valid for any beam whose J z is well defined. For the specific case of an 
incident paraxial LG p z beam, the integral in eq. © can be solved analytically and the following 
expression is obtained: 

Cjip = [U +p)Kj -pY-U + i+pV-U - i-pW 2 

(-iy M l+2s (f qtl ) (4) 



E 



(j - l - p - s)\(j +p - s)\(s + l)\s\ 2 l + 2s 



where the number M a (f ql ) = J °° k r dk r k^f q ^{k r ) is the momentum of order a of the function f qi . 
The expression for f q ,i(k r ) when the incident beam is a LGg^ is: 



fM = ( 27r( ^ g)! ) V2 (WW) L^{wlkl/2) exp {-wlkl/A) (5) 

with wq being the beam width of the LG mode in real space and L l q (x) the associated Laguerre 
polynomials. Finally, for completeness sake, we give the expression of M a (f qt i) when the function 
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FIG. 1. \Cji p \ 2 for different LG P /. a) LGo,o b) LGo,2 c) LGo,4 d) LGo,6- The relation between the width (w) 
of the beam and the wavelength is wj A = 4. 

f Qt i is given by the expression © and we also particularize the result for q = 0: 
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where 2-^1(0, b,c: z) is the hypergeometric function. The problem, then, is solved. Using eqs. ([3]), 
© and © the decomposition of the incident beam is found and then the scattered and the interior 
fields are obtained by applying Maxwell boundary conditions. As a result, we obtain: 



j=\i+p\ 



.1 



E ^(2j + l) 1/2 ^ p [c,A^ p) + zHA 
j=\i+ P \ 



(e) 



(8) 



where Cj| p is given by eq. © and its shape is depicted in Fig.l for different cases. Also, 
{a n , b n , c n , d n } are the classical Mie coefficients defined in lil8L Chap. 4]. Note that the two 
fields have the same formal expression as the incident field, except for the Mie coefficients and the 
radial dependence of the mulipolar fields, which is a Bessel function for the incident and interior 
fields and a Hankel function for the scattered one. These coefficients only depend on the radius 
of the particle (R), and the relative permeability (/x r ) and permittivity (e r ) of the sphere with re- 
spect to the surrounding medium. Next, the analytical expressions for aj and bj are provided when 



fi r = 1: 



n r ipj (n r x)ip'j (x) — ipj (x)if>j (n r x) ipj (n r x)ipj (x) — n r ipj {x)if>'j (n r x) 



a . = ' ' T J K ' TJK ,T J K ' ' jj . = TJX ' ' T ^ K ' ' TJK illl ' ' (9) 

n r ijjj{nrx)^{x) - ^(xWjforx) 1 il)j{n r x)Q{x) — n r ^j{x)ijj'j{n r x) 

where ipj and ^ are the Riccati-Bessel functions of order j and x = 2nr /A is the so-called size 

parameter, and where n r = y/fi r e r . 



CROSS SECTIONS 

Once the fields have been obtained, the scattering efficiency associated with them can be cal- 
culated. This efficiency is related to the power of a beam in the far field. Its definition, for the case 
of a plane wave excitation, is given by the flux of the Poynting vector across a spherical surface 
divided by the product of the incident irradiance and the particle's cross-sectional area projected 



onto a plane perpendicular to the incident beam II18L Chap. 3]: 



Qme = j-^p j S sm ■ dSl (10) 

with S sca = Re {(l/2)E sca x H sca *} the Poynting vector, dVL an element of surface perpendicular 
to the surface, and Ii the incident irradiance. This expression has a major drawback when it is 
applied to incident fields that vary point to point. Hence, we have used another normalization 
factor to compute the scattering efficiency. We normalized the expressions over the integral of 
|E*| 2 across a transverse surface: 

Q sca = 1 f S sca ■ dtt (11) 



To perform this calculation, we must derive the electric and magnetic fields from the vector poten- 
tial in a Coulomb gauge, that is: 

E = ikA (12) 
H = VxA (13) 

Then, the integral (fTTj) is performed paying special attention to the orthonormality relations of the 
multipolar fields, which are: 

" AHxAW').rdn = (14) 



/ 



AgxAgH-r ,/(> = () (15) 

dQ A^Ag?* = F (r) 5 jf 8 mm ,5 xx , (16) 



where F(r) is 



y 2 

= < jio + ii 2 +(j + i)i^-ii 



for x 



for a; = e 



(17) 



2j + l 

and £ 3 - is a Bessel function for the incident and interior field and a Hankel function of the first kind 
for the scattered field. 

Finally, the scattering efficiency for the incident beam in eq. © is: 



Q" 



j=\i+p\ 



27T (2j 



A- 2 



|CjZp| 2 (l a i| 2 



+ N 2 ) 



(18) 



with Cj7 P given by eq. © and where the following relation holds + l)\Cji p \ 2 = 2 if the 

incident field is normalized to 1 |l7|1 . In contradiction, the result from the Mie theory [7] is: 



Q 



sea 
Mie 



(19) 



Expressions (1181) and (1191) give rise to very different scattered cross sections. We show these 
differences in Fig. 2, 3 and 4. First, expression (fT9l) is used to compute the Mie scattering cross 
section. The result of this computation is depicted in Fig. 2. It can be seen that the scattering cross 
section behaviour is dominated by the low order modes and that the ripples are caused by the high 
order electromagnetic modes. Next, expression (fT8l) is used to plot Fig. 3 and Fig. 4. In Fig. 3 the 




FIG. 2. Q sca and some of the Mie coefficients for the conventional Mie Theory, when m = 1.33 + HO 8 



and x £ [50, 52] (Ref. UM, Fig. 11.7) 



scattering cross section is plotted for four different LG g; ; modes. It can be observed that the angular 
momentum of the incident beam plays a crucial role. Indeed, the ripple structure is increasingly 
enhanced as the angular momentum of the LG beam in consideration is increased. This effect 



7 



is specially clear in Fig. 3d where some of the ripples present in Fig. 3a are largely enhanced. 
Finally, Fig. 4 shows the enhancement of a single ripple. Although presented in a different way, 
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FIG. 3. Q sca as a function of the wavelength for different incident LG beams, a) LGo,o b) LGo,2 c) LGo,4 
d) LGo,6- The other parameters are kept equal to n r = 1.5, p = 1, r = 1.3 \xm. The decomposition of each 
of the LG beams used to computed these cross sections is shown in Fig. 1. 
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FIG. 4. Q sca for a very short range of wavelegths. a) LGo,o b) LGo,2 c) LGo,4 d) LGo,6- Equally to Fig. 3, 
the other variables are kept equal ton r = 1.5, p = l,r = 1.3 \im. 

the same effect presented in Fig. 3 is observed again. However, in this case, the enhancement of 
the ripple structure can be determined quantitatively. We have compared the maximum value of 
the cross section which peaks at 560.3 nm and the minimum value. This minimum value has been 
computed as the average of two values at 0.7nm apart in wavelength. Then, we have computed the 
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Fig. 4 LG ,o LG 0i 2 LG ,4 LG 0) 6 
7 (%) 1.455 6.916 17.30 113.6 

TABLE I. Enhancement of the ripple structure using different LG beams. 



ratio 7, which is the relative quotient between the subtraction of the maximum and the minimum, 
and the minimum. 

Qmax Qmin zor>\ 

7 = 7^ (20) 

where Q max = QIqq 3 and Q m in = (Qtei o + Q559 e)/^- The results are presented in Table 1. It can 
be seen that the background is highly reduced. When a LG ,o is used to excite the sphere, the peak 
only represents 1% of the background. Nonetheless, when a LG ,6 is used, the peak contribution 
is larger than the background's one. That is, the background has been reduced more than a half. 



CONCLUSIONS 



We have shown that the LG beams can be used to enhance the ripple sctructure in dielec- 
tric spheres. The conservation of the AM of light implies that the first multipolar modes do not 
contribute to the scattering of the particle and therefore the ripple structure can be turned into res- 
onances. This is a great improvement towards the control of scattered fields. This new effect could 
be applied in a large variety of fields such as spectroscopy, cytometry or dark-field microscopy to 
mention a few. 
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